On Non-Gaussianity in the Curvaton Scenario 
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Since a positive future detection of non-linearity in the cosmic microwave background anisotropy 
pattern might allow to descriminate among different mechanisms giving rise to cosmological adia- 
batic perturbations, we study the evolution of the second-order cosmological curvature perturbation 
on super-horizon scales in the curvaton scenario. We provide the exact expression for the non- 
Gaussianity in the primordial perturbations including gravitational second-order corrections which 
are particularly relevant in the case in which the curvaton dominates the energy density before it 
decays. As a byproduct, we show that in the standard scenario where cosmological curvature per- 
turbations are induced by the inflaton field, the second-order curvature perturbation is conserved 
even during the reheating stage after inflation. 
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I. INTRODUCTION 

One of the basic ideas of modern cosmology is that 
there was an epoch early in the history of the universe 
when potential, or vacuum, energy associated to a scalar 
field, the inflaton, dominated other forms of energy den- 
sity such as matter or radiation. During such a vacuum- 
dominated era the scale factor grew exponentially (or 
nearly exponentially) in time. During this phase, dubbed 
inflation [1,2], a small, smooth spatial region of size less 
than the Hubble radius could grow so large as to easily 
encompass the comoving volume of the entire presently 
observable universe. If the universe underwent such a 
period of rapid expansion, one can understand why the 
observed universe is so homogeneous and isotropic to high 
accuracy. 

Inflation has also become the dominant paradigm for 
understanding the initial conditions for structure for- 
mation and for Cosmic Microwave Background (CMB) 
anisotropy. In the inflationary picture, primordial den- 
sity and gravity-wave fluctuations are created from quan- 
tum fluctuations "redshifted" out of the horizon during 
an early period of superluminal expansion of the uni- 
verse, where they are "frozen" [3-7] . Perturbations at the 
surface of last scattering are observable as temperature 
anisotropy in the CMB, which was first detected by the 
Cosmic Background Explorer (COBE) satellite [8-10]. 
The last and most impressive confirmation of the infla- 
tionary paradigm has been recently provided by the data 
of the Wilkinson Microwave Anistropy Probe (WMAP) 
mission which has marked the beginning of the preci- 
sion era of the CMB measurements in space [11]. The 
WMAP collaboration has produced a full-sky map of the 
angular variations of the CMB, with unprecedented ac- 
curacy. WMAP data confirm the inflationary mechanism 
as responsible for the generation of curvature (adiabatic) 
super-horizon fluctuations. 

Despite the simplicity of the inflationary paradigm, the 



mechanism by which cosmological adiabatic perturba- 
tions are generated is not yet fully established. In the 
standard picture, the observed density perturbations are 
due to fluctuations of the inflaton field itself. When in- 
flation ends, the inflaton oscillates about the minimum of 
its potential and decays, thereby reheating the universe. 
As a result of the fluctuations each region of the universe 
goes through the same history but at slightly different 
times. The final temperature anisotropics are caused by 
the fact that inflation lasts different amounts of time in 
different regions of the universe leading to adiabatic per- 
turbations. Under this hypothesis, the WMAP dataset 
already allows to extract the parameters relevant for dis- 
tinguishing among single- field inflation models [12]. 

An alternative to the standard scenario is represented 
by the curvaton mechanism [13-16] where the final curva- 
ture perturbations are produced from an initial isocurva- 
ture perturbation associated to the quantum fluctuations 
of a light scalar field (other than the inflaton), the curva- 
ton, whose energy density is negligible during inflation. 
The curvaton isocurvature perturbations are transformed 
into adiabatic ones when the curvaton decays into radi- 
ation much after the end of inflation*. Contrary to the 
standard picture, the curvaton mechanism exploits the 
fact that the total curvature perturbation (on uniform 
density hypersurfaces) Q can change on arbitrarily large 
scales due to a non-adiabatic pressure perturbation which 
may be present in a multi-fluid system [18-22]. While the 
entropy perturbations evolve independently of the cur- 



* Recently, another mechanism for the generation of cosmo- 
logical perturbations has been proposed [17]. It acts during 
the reheating stage after inflation if super-horizon spatial fluc- 
tuations in the decay rate of the inflaton field are induced 
during inflation, causing adiabatic perturbations in the final 
reheating temperature in different regions of the universe. 
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vature perturbation on large scales, the evolution of the 
large-scale curvature is sourced by entropy perturbations. 

Fortunately the standard and the curvaton scenarios 
have different observational signatures. The curvaton 
scenario allows to generate the observed level of den- 
sity perturbations with a much lower scale of inflation 
and thus generically predicts a smaller level of gravita- 
tional waves. More interestingly, density perturbations 
generated through the curvaton scenario could be highly 
non-Gaussian and the level of non-Gaussianity in the pri- 
mordial perturbations, which is usually parametrized by 
a dimensionless non- linear parameter /nl, depends upon 
an unknown parameter r indicating the fraction of energy 
density contributed by the curvaton field at the epoch of 
its decay. For tiny values of r, it has been estimated that 
the non-Gaussianity can be large enough to be detectable 
by present CMB experiments [14,16]; the current WMAP 
[23] bound on non-Gaussianity, |/nl| < 10 2 , already re- 
quires r to be larger than about 10~ 2 . 

This curvaton prediction has to be contrasted to what 
predicted within the traditional one-single field model of 
inflation where the initially tiny non-linearity in the cos- 
mological perturbations generated during the inflation- 
ary epoch [24,25] gets enhanced in the post-inflationary 
stages giving rise to a well-defined prediction for the non- 
linearity in the gravitational potentials [26] . 

Since a positive future detection of non-linearity in 
the CMB anisotropy pattern might allow to descriminate 
among the mechanisms by which cosmological adiabatic 
perturbations are generated, it is clear that the precise 
determination of the non-Gaussianity predicted by the 
curvaton mechanism is of primary interest. 

The goal of this paper is to provide an exact expres- 
sion for the non-linear parameter /nl within the curvaton 
scenario including second-order corrections from gravity 
which are particularly relevant in the case in which the 
curvaton dominates the energy density before it decays. 
We perform a fully relativistic analysis of the dynam- 
ics of second-order perturbations taking advantage of the 
second-order gauge-invariant curvature perturbation in- 
troduced in Refs. [27,28] (see also [24,29,30]) and showing 
how it evolves on arbitrarily large scales in the presence 
of two fluids, matter (the curvaton) and radiation. Our 
results generalize the estimates given in Refs. [14,16] and 
confirm their findings in the limit r<l 

The paper is organized as follows. In Section II we 
briefly summarize the properties of the curvaton sce- 
nario and how the primordial curvature perturbations 
are created at first-order. In section III we compute the 
second-order curvature perturbation from the curvaton 
fluctuations and determine the exact expression for the 
non-linear parameter /nl as a function of the unknown 
parameter r. In section IV we show that our findings 
can be easily generalized to the standard scenario where 
adiabatic perturbations are provided by the same field 
driving inflation and prove that the second-order curva- 
ture perturbation is conserved even during the reheating 



stage after inflation when the inflaton field decays to give 
birth to the standard radiation phase. Finally, Section V 
contains our conclusions. 



II. GENERATING THE CURVATURE 
PERTURBATION AT LINEAR ORDER 

During inflation the curvaton field a is supposed to 
give a negligible contribution to the energy density and 
to be an almost free scalar field, with a small effective 
mass m 2 = \d 2 V/da 2 \ <C iff [14,16], where if/ = a/a is 
the Hubble rate during inflation, a is the scale factor and 
a dot denotes derivative with respect to cosmic time. 

The unperturbed curvaton field satisfies the equation 
of motion 

a" + 2Ha' + a 2 ^=0, (1) 

OCT 

where a prime denotes diffcrcntation with respect to the 
conformal time dr — dt/a and H = a' /a is the Hub- 
ble parameter in conformal time. It is also usually as- 
sumed that the curvaton field is very weakly coupled to 
the scalar fields driving inflation and that the curvature 
perturbation from the inflaton fluctuations is negligible 
[14,16]. Thus, if we expand the curvaton field up to first- 
order in the perturbations around the homogeneous back- 
ground as <t(t, x) — a(r) + S^ a, the linear perturbations 
satisfy on large scales 

5^a" + 2H5^a' + a 2 ^5^a = 0. (2) 

As a result on superhorizon scales its fluctuations 5a will 
be Gaussian distributed and with a nearly scale- invariant 
spectrum given by 

f^, (3) 

where the subscript * denotes the epoch of horizon exit 
k = aH. Once inflation is over the inflaton energy den- 
sity will be converted to radiation (7) and the curvaton 
field will remain approximately constant until if 2 ~ m 2 . 
At this epoch the curvaton field begins to oscillate around 
the minimum of its potential which can be safely approx- 
imated to be quadratic V w ^m 2 a 2 . During this stage 
the energy density of the curvaton field just scales as non- 
rclativistic matter p a cx a~ 3 . The energy density in the 
oscillating field is 

Pa(r,x) w m 2 cr 2 (r,x) , (4) 

and it can be expanded into a homogeneous background 
Pct(t) and a first-order perturbation 8^p a as 

p CT (r,x) = p a {r) + (5 (1) p<t(t,x) = m 2 cr + 2m 2 a <5 (1) ct . 

(5) 
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As it follows from Eqs. (1) and (2) for a quadratic poten- 
tial the ratio S^'cr/cr remains constant and the resulting 
relative energy density perturbation is 



= 2 fs^'o 



a 



(6) 



where the * stands for the value at horizon crossing. 

Such perturbations in the energy density of the cur- 
vaton field produce in fact a primordial density pertur- 
bation well after the end of inflation. The primordial 
adiabatic density perturbation is associated with a per- 
turbation in the spatial curvature ip an d it is usually 
characterized in a gauge-invariant manner by the cur- 
vature perturbation £ on hypersurfaces of uniform total 
density p. At linear order the quantity ( is given by the 
gauge- invariant formula [31] 



P 1 



(7) 



and on large scales it obeys the equation of motion [31,21] 



C (D' = _ 



H 



P + P 



nad 



where c5«P nad = S^P - c 2 s 5^p is the non-adiabatic 
pressure perturbation, S^P being the pressure pertur- 
bation and c 2 s = P' I p' the adiabatic sound speed. In the 
curvaton scenario the curvature perturbation is gener- 
ated well after the end of inflation during the oscillations 
of the curvaton field because the pressure of the mixture 
of matter (curvaton) and radiation produced by the in- 
flaton decay is not adiabatic. A convenient way to study 
this mechanism is to consider the curvature perturba- 
tions Q associated with each individual energy density 
components, which to linear order are defined as [21] 



(9) 



Therefore, during the oscillations of the curvaton field, 
the total curvature perturbation in Eq. (7) can be written 
as a weighted sum of the single curvature perturbations 
[21,16] 



C (1) = (1 - + f& 
where the quantity 

3pa 



f = 



4p 7 + 3p CT 



(10) 



(11) 



defines the relative contribution of the curvaton field to 
the total curvature perturbation. From now on we shall 
work under the approximation of sudden decay of the cur- 
vaton field. Under this approximation the curvaton and 
the radiation components p a and p 7 satisfy separately 
the energy conservation equations 



p 7 = -m Pl , 

Pa = , 



(12) 



and the curvature perturbations Q remains constant 
on supcrhorizon scales until the decay of the curvaton. 
Therefore from Eq. (10) it follows that the first-oder cur- 
vature pertubation evolves on large scales as 

C (1) ' - /'(C< 1} - C 7 X) ) = Hf(l /)(#> - C 7 1} ) , (13) 

and by comparison with Eq. (8) one obtains the expres- 
sion for the non-adiabatic pressure perturbation at first 
order [14,16] 



^P„ad=P.(l-/)(d 1) -d 1) ). 



(14) 



Since in the curvaton scenario it is supposed that the 
curvature perturbation in the radiation produced at the 
end of inflation is negligible 



7 4 p 7 



(8) Similarly the value of ^ is fixed by the fluctuations of 



0. 



(15) 



the curvaton during inflation 



" Pa 



(16) 



where I stands for the value of the fluctuations during 
inflation. From Eq. (10) the total curvature perturbation 
during the curvaton oscillations is given by 



C (1) = /#> • 



(17) 



As it is clear from Eq. (17) initially, when the curvaton 
energy density is subdominant, the density perturbation 
in the curvaton field gives a negligible contribution to 
the total curvature perturbation, thus corresponding to 
an isocurvature (or entropy) perturbation. On the other 
hand during the oscillations p a oc a~ 3 increases with re- 
spect to the energy density of radiation p 7 oc a~ 4 , and 
the perturbations in the curvaton field are then converted 
into the curvature perturbation. Well after the decay of 
the curvaton, during the conventional radiation and mat- 
ter dominated eras, the total curvature perturbation will 
remain constant on superhorizon scales at a value which, 
in the sudden decay approximation, is fixed by Eq. (17) 
at the epoch of curvaton decay 



C (1) = .feC 



(1) 



(18) 



where D stands for the epoch of the curvaton decay. 

Going beyond the sudden decay approximation it is 
possible to introduce a transfer parameter r defined as 
[16,22] 



C (D = rC (i 



(19) 
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where is evaluated well after the epoch of the cur- 
vaton decay and is evaluated well before this epoch. 
The numerical study of the coupled perturbation equa- 
tions has been performed in Ref. [22] showing that the 
sudden decay approximation is exact when the curvaton 
dominates the energy density before it decays (r = 1), 
while in the opposite case 



Pa 

P / D 



III. SECOND-ORDER CURVATURE 
PERTURBATION FROM THE CURVATON 
FLUCTUATIONS 

Here we generalize to second-order in the density per- 
turbations the results of the previous section. 

As it has been shown in Ref. [28] it is possible to define 
the second-order curvature perturbation on uniform total 
density hypersurfaces by the quantity (up to a gradient 
term) 



c 



(2) 



+ 2H 



p> pi 



2 i_£(yi)' + 2wv< (1) ) 



8W, 



H^-r - H' - 2H 2 
P 



(21) 



where the curvature perturbation ip has been expanded 
up to second-order as ip = ip^ + ^ip^ and 8^p cor- 
responds to the second-order perturbation in the total 
energy density around the homogeneous background p(r) 

p{r,x) = P (t) + 5p(T,x) 

= p(T)+<5«p(T,x) + ^ 2 )p(r,x). (22) 

The quantity £( 2 ) is gauge-invariant and, as its first-order 
counterpart defined in Eq. (7), it is sourccd on super- 
horizon scales by a second-order non-adiabatic pressure 
perturbation [28]. 

In Ref. [26] the conserved quantity has been used 
in the standard scenario where the generation of cosmo- 
logical perturbations is induced by fluctuations of the 
inflaton field (and there is no curvaton) in order to fol- 
low the evolution on large scales of the primordial non- 
linearity in the cosmological perturbations from a period 
inflation to the matter dominate era. In the present sce- 
nario the conversion of the curvaton isocurvature pertur- 
bations into a final curvature perturbation at the epoch 
of the curvaton decay can be followed through the sum 
(10) of the individual curvature perturbations weighted 
by the ratio / of Eq.(ll). 



Let us now extend such a result at second order in 
the perturbations. As we shall sec in Sec. Ill A this re- 
sult will enable us to compute in an exact way the level 
of non-Gaussianity produced by the non-linearity of the 
perturbations in the curvaton energy density. 

Since the quantities (j 1 ^ and C^" 1 are gauge-invariant, 
we choose to work in the spatially flat gauge ip = if not 
otherwise specified. Note that from Eqs. (6) and (16) the 
(i) 



(20) value of Q is thus given by 



u (1 V 

3 Pa 



2 6^a 

3 (T 



5W* 



(23) 



where we have used the fact that £^ ( or equivalently 
S^a/a) remains constant, while from Eq. (16) in the 
spatially flat gauge 



AD _ 1 6 {1) Pj 
^ 4 Pl 



(24) 



During the oscillations of the curvaton field the first-order 
energy conservation equations in the spatially flat gauge 
ip = yield on large scales 1. 



= 5^ P ' a + a(V 7 = -ms^pv - 4^«p 7 . 

and hence using Eqs. (12), (23), (24) and (25) 



(25) 



H 



P' 
5Wp 



3/#>+ 4(1 -/)#>, 



>(i) 

^7 



(26) 



We can thus rewrite the total second-order curvature per- 
turbation £( 2 ) as 



c (2) 



-H 



5^h 



Ka^ + (1 - 



/ 2 d 1) + (i-/)(2 + /)d 1) 



(27) 



In a similar manner to the linear order, let us introduce 

(2) 

now the curvature perturbations Q at second order for 
each individual component. Such quantities will be given 
by the same formula as Eq. (21) relatively to each energy 
density pi 



^Here and in the following we neglect gradient terms which, 
upon integration over time, may give rise to non-local op- 
erators which are not necessarily suppressed on large-scales 
being of the form V~ 2 [V(-)V(-)] or V~ 2 [(-)V 2 (-)]- In this 
paper we will focus on the momentum-independent part of 
the non-linear parameter /nl- 
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42) 



6^ pi 



5 W 



Pi 



2 d —Pl ^(i)' + 2WV (1) ) 
H' -2H 2 ) . (28) 



Using the same procedure described above it follows 
that in the spatially flat gauge 



A2)_ 

W 3 Pa 



r( 2 ) = ± 



Pi 



- (tf ') 2 • 



(29) 
(30) 



Such quantities are gauge-invariant and, in the sudden 
decay approximation they are separately conserved until 
the curvaton decay. Using Eqs. (29) and (30) to express 
the second-order perturbation in the total energy density 
S^p = S^pa+S^p-y, and after some algebra, one finally 
obtains the following expression for the total curvature 
perturbation 



C (2) = ft™ + (i - /Kf 



(31) 



Eq. (31) is one of our main results. It generalizes to 
second-order in the perturbations the weighted sum of 
Eq. (10). In particular notice that in the limit where 
one of the two fluid is completely subdominant (/ — > 

(2) 

or / — > 1) the corresponding curvature perturbation Q 
turns out to coincide with the total one (^ 2 ) . 

Under the sudden decay approximation of the curva- 
ton field the individual curvature perturbations are sep- 
arately conserved on large scales, and thus from Eq. (31) 
it follows that C*- 2 -* evolves according to the equation 

c (2) ' = /' (6 2) cf ) + /'(i - 3/ 2 ) (c« - c< 1} ) 2 • 

(32) 

Note that Eq. (32) can be rewritten as [28] 

n 



£(2)' 



p + P 

2 



-(5( 2 )P 



< 5( 1 )p nad -2(p + P)C (1) lc (1) ', (33) 



P + P 

with J^Pnad given by Eq. (14) and 

= p„{\- /) [ (cf - 6 2) ) + if + 6/ - 1) 

x^-C^' + ^^-Cf)], (34) 

is the gauge-invariant non-adiabatic pressure perturba- 
tion on uniform density hypersurfaces on large scales 



which can be checked to coincide with the generic ex- 
pression provided in Ref. [28] 



§(2)p = S (2)p_^ S (2) p + pl 
P 



pf P' 



P p 



~7 



p' 

(35) 



The second-order curvature perturbation in the standard 
radiation or matter eras will remain constant on super- 
horizon scales and, in the sudden decay approximation, 
it is thus given by the quantity in Eq. (31) evaluated at 
the epoch of the curvaton decay 



c (2) = /z4 2) + Mi-/!>) (c«y 



(36) 



where we have used the curvaton hypothesis that the 
curvature perturbation in the radiation produced at the 
end of inflation is negligible so that (£p ~ and ffl w 
0. The curvature perturbation <^ is given by Eq. (23), 

(2) 

while (a in Eq. (29) is obtained by expanding the energy 
density of the curvaton field, Eq. (4) , up to second order 
in the curvaton fluctuations 

A,(x,i) - pAr) + 5^p a {r,x i ) + l -5^ Pa (T,x l ) 

= mla + 2mlaS^a + ml (^V)* . (37) 



It follows that 



5^ Pa = 1 (6V> Po 

Pn 2 



-§(<?>)'• 



(38) 



where we have used Eq. (23), and hence from Eq. (29) 
we obtain 



(39) 



(2) 

where we have emphasized that also Q is a conserved 
quantity whose value is determined by the curvaton fluc- 
tuations during inflation. Plugging Eq. (39) into Eq. (36) 
the curvature perturbation during the standard radiation 
or matter dominated eras turns out to be 



A. Non-Gaussianity of the curvaton perturbations 

Let us now now focus on the calculation of the non- 
linear parameter /nl which is usually adopted to char- 
acterized the level of non-Gaussianity of the Bardeen po- 
tential [32]. In order to compute / NL in the curvaton 



•5 



scenario, we switch from the spatially flat gauge ip = 
to the longitudinal or Poisson gauge [33] . Such a proce- 

(2) 

dure is possible since the curvature perturbations Q are 
gauge- invariant quantities. In particular this is evident 
form the expression found in Eq. (40). During the matter 
dominated era from Eq. (21) it turns out that [26] 



C (2) = _^(2) + 



+ 



3 p 



1 6 - r> 
3 P 
1 S^p 



where in the last step we have used that on large scales 
S^p/p = -2-0W in the Poisson gauge [26]. Eq. (41) 
combined with Eq. (40), which gives the constant value 
on large scales of the curvature perturbation C^ 2 -* during 
the matter dominated era, yields 



1(5( 2 ) 



3 p 



75 

20 - — + 25 f D 



(42) 



where we have used JdC^ = — ^ip^ from Eq. (18) and 
the usual linear relation between the curvature pertur- 
bation and the Bardeen potential C^ 1 ) = — |V ,< - 1 - ) during 
the matter dominated era. Since on large scales (from 
the second-order (0 — 0)- and (i — j)-components of Ein- 
stein equations, see Eqs. (A. 39) and (A. 42-43) in [24]) 
the following relations hold during the matter-dominated 
phase 



9 2 p 



+4(> i )y 



- 10 V" 2 (&dj , (43) 

we conclude that 



(2) = 



10 5 

y + 3 /D 



5 



(0(ki)0(k 2 )0(k3)) = {2irf (ki + k 2 + k 3 ) 
2/* L (ki,k 2 ) VtikJV+fa) + cyclic 



(46) 



where P,p(k) is the power-spectrum of the gravitational 
potential. From Eq. (44) we read the the non-linearity 
parameter 



NL 
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5 

4r 



5(ki,k 2 ) 



(41) where 



(44) 



The total curvature perturbation will then have a non- 
Gaussian (x 2 )_ com p 0ne nt. The lapse function <j> = (j)^ + 
can be expressed in momentum space as 

4>(k) = <^ (1) (k) + ^3 J d 3 h! d 3 k 2 <5 (3) (ki + k 2 - k) 
x/^ L (k 1 ,k 2 )^ 1 )(k 1 )^ 1 )(k 2 ), (45) 

where we have defined an effective "momentum- 
dependent" non-linearity parameter /^ L . Here the linear 
lapse function (jP^ — ip^ is a Gaussian random field. 
The gravitational potential bispectrum reads 



ff(ki,k 2 ) 



ki ■ k 2 

fc 2 



1 + 3 



kl k; 

k 2 



(47) 



(48) 



with k = ki + k 2 and we have replaced fo with r to 
go beyond the sudden approximation. Notice that in the 
final bispectrum expression, the diverging terms arising 
from the infrared behaviour of /^ L (k 1? k 2 ) are automat- 
ically regularized once the monopole term is subtracted 
from the definition of <f> (by requiring that (<f>)=0). 

As far as the momentum-independent part is con- 
cerned we note that in the limit r <C 1 we obtain 
/nl = — T" wn ich reproduces the estimate provided in 
[14,16], while, in the limit r ~ 1, we obtain = | * for 
r ~ 1. These values have to be compared to the value 
/nl = ~~ \ [26] obtained for perturbations whose wave- 
lengths re-enter the horizon during the matter-dominated 
phase for the standard scenario in which curvature per- 
turbations are induced by fluctuations of the inflaton 
field. We conclude that if r <C 1 the non-Gaussianity 
in the curvaton scenario is larger than the one pre- 
dicted in the standard scenario. Finally we point out 
that additional non-Gaussianity will be generated after 
horizon-crossing, due to known Newtonian and relativis- 
tic second-order contributions which are relevant on sub- 
horizon scales, such as the Rees-Sciama effect [34], whose 
detailed analysis has been given in Refs. [35]. It is impor- 
tant to consider also these effects when making a com- 
parison with observations. 



IV. A COMMENT ON THE EVOLUTION OF 
THE CURVATURE PERTURBATION DURING 
THE REHEATING PHASE IN THE STANDARD 
SCENARIO 

Note that Eq. (31) is indeed valid in the general frame- 
work of an oscillating scalar field and a radiation fluid, 
the curvaton scenario being only a particular case. Thus 
in this section we shall indicate the generic scalar field 



^Notice that the formula (36) in [16] for the estimate of the 
non-linear parameter contains a sign misprint and should read 

/nl - -£> S ivin S /nl 



- 1 for r ■ 



1. 
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by (p instead of a. From Eq. (31) it is possible to derive 
the following equation of motion on large scales for the 
second-order curvature perturbation C}^ 



the curvature perturbations are generated by spatial vari- 
ations of the inflaton decay rate [17]. 



c (2) ' = r (4 2) - cf ) + /'(i - 3.n (#> - #>) , 

(49) 

where we have used the fact that, in the approximation 
of sudden decay of the scalar field ip, the individual cur- 
vature perturbations at first and second-order are sepa- 
rately conserved. Using Eqs. (10) and (31) it is possible 

to rewrite C (2) ' in terms only of C (2) , c£° and C (1) , 
as 

c (2) ' = -Hf (c (2) - 4 2) ) + Hf{\ + 2f) (c« - c«) 2 . 

(50) 

Here we want to make a simple but important observa- 
tion. Besides the curvaton scenario, the most interesting 
case where there is an oscillating scalar field and a ra- 
diation fluid is just the phase of reheating following a 
period of inflation in the standard scenario for the gen- 
eration of cosmological perturbations on large scales. In 
such a situation the oscillating scalar field is just the in- 
flaton field ip whose fluctuations induce curvature per- 
turbations. Therefore it is possible to see in a straight- 
forward way that during the reheating phase, when the 
inflaton field finally decays into radiation, a solution of 
Eq. (50) is the one corresponding to a total curvature 
perturbation which is indeed fixed by the inflaton cur- 
vature perturbation during inflation £W = = CiV> 

V. CONCLUSIONS 

In this paper we have determined the evolution on 
large scales of the second-order curvature perturbation 
within the curvaton scenario where two fluids are present, 
the curvaton and radiation. We have computed the the 
non-linear parameter /nl measuring the level of non- 
Gaussianity in the primordial cosmological perturbations 
and provide its exact expression as a function of the pa- 
rameter r. Our findings are particularly interesting if 
one wishes to extract from a positive future detection of 
non-linearity in the CMB anisotropy pattern a way to de- 
scriminate among the mechanisms by which cosmological 
adiabatic perturbations are generated. It would be inter- 
esting to extend our results to those models which can ac- 
comodate for a primordial value of /nl larger than unity. 
This is the case, for instance, of a large class of multi- 
field inflation models which leads to cither non-Gaussian 
isocurvature perturbations [36] or cross-correlated non- 
Gaussian adiabatic and isocurvature modes [37] and the 
so-called "inhomogeneous reheating" mechanism where 
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